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Bateman's principles explain sex roles and sexual dimorphism through 
sex-specific variance in mating success, reproductive success and their 
relationships within sexes (Bateman gradients). Empirical tests of these prin- 
ciples, however, have come under intense scrutiny. Here, we experimentally 
show that in replicate groups of red junglefowl, Gallus gallus, mating and 
reproductive successes were more variable in males than in females, resulting 
in a steeper male Bateman gradient, consistent with Bateman's principles. 
However, we use novel quantitative techniques to reveal that current methods 
typically overestimate Bateman's principles because they (i) infer mating suc- 
cess indirectly from offspring parentage, and thus miss matings that fail to 
result in fertilization, and (ii) measure Bateman gradients through the univari- 
ate regression of reproductive over mating success, without considering 
the substantial influence of other components of male reproductive success, 
namely female fecundity and paternity share. We also find a significant 
female Bateman gradient but show that this likely emerges as spurious con- 
sequences of male preference for fecund females, emphasizing the need 
for experimental approaches to establish the causal relationship between 
reproductive and mating success. While providing qualitative support for 
Bateman's principles, our study demonstrates how current approaches can 
generate a misleading view of sex differences and roles. 



Author for correspondence: 

Julie M. Collet 

e-mail: j.collet@uq.edu.au 



Electronic supplementary material is available 
at http://dx.doi.org/10.1098/rspb.2013.2973 or 
via http://rspb.royalsocietypublishing.org. 



1. Introduction 

In a pioneering study published in 1948, Bateman [1] extrapolated from experimen- 
tal results in Drosophila melanogaster to propose that intrasexual selection is normally 
more intense in males because t3^ically (i) compared with females, males 
have higher variance in number of mates (i.e. mating success); (ii) males also have 
higher individual variation in the number of offspring produced (i.e. reproductive 
success) than females and (iii) the slope of the relationship between mating and 
reproductive success (Bateman gradient) is steeper in males than in females. 
These observations became known as Bateman's principles [2,3]; their formalization 
through selection analysis [4,5] marked the advent of modem sexual selection 
theory [5-7]. Bateman's principles provide a conceptual explanation for Darwin's 
observations that sexual selection is typically more intense in males and that 
males are often 'eager' to mate whereas females are 'coy' [8]. This Darwin- Bateman 
paradigm represents the foundation of our understanding of the evolutionary 
ecology of sex-specific selection, sex roles and sexual dimorphism [3,9]. 

A number of empirical studies have provided qualitative support for 
Bateman's principles and their link with sex roles [3,10-17]. Recent work, how- 
ever, has highlighted a number of problems with both the measurement and 
interpretation of Bateman's principles [18-21], calling into question the magnitude 
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of the estimated sex differences and their eco-evolutionary 
significance for four main reasons. First, Bateman gradients 
are measured by inferring mating success from offspring 
parentage [5,19,22,23], based on the assumption that an individ- 
ual sires offspring with each of its mates. However, this 
assumption is problematic especially when polyandry gener- 
ates sperm competition, which can result in males gaining no 
paternity despite behaviourally successful matings. Inferring 
mating success from offspring parentage in the absence of inde- 
pendent data on mating behaviour can thus introduce a source 
of error in estimates of male Bateman gradients. In most 
studies, measures of mating success based on parentage also 
capture, implicitly or explicitly, components of offspring survi- 
val (e.g. [24], reviewed in reference [25]), which can further bias 
estimates of Bateman gradients [26,27]. However, the magni- 
tude of this potential source of error remains untested. 
Moreover, because — all else being equal — the probability of a 
male fertilizing at least one egg is a function of clutch size, 
the bias introduced by inferring male mating success from 
parentage is expected to be inversely proportional to average 
clutch size [22]. A better alternative to inferring mating success 
from offspring parentage is to observe mating behaviour 
directly and record mating rates for each individual [14]. How- 
ever, the use of mating rates often fails to distinguish between 
an individual mating multiple times with the same partner, or 
mating once with multiple partners [28]. This distinction is cru- 
cial, because remating with the same partner and mating with 
multiple partners may influence the reproductive success of an 
individual in drastically different ways. In males, for example, 
remating with the same female would influence a male's 
chances of fertilizing her eggs, whereas mating with multiple 
females would directly contribute to his mating success. One 
would therefore need to measure the actual number of mates 
for each individual, controlling for the potentially confounding 
effect of number of matings obtained with each mate. This 
requires intensive observations of individually tagged males 
and females [17,29]. 

Second, measuring the Bateman gradient as the slope of the 
linear least-square regression of reproductive success over 
mating success may be simplistic, because even with unlimited 
availability of mating opportunities the reproductive success 
of an individual is, at some point, limited by intrinsic repro- 
ductive costs such as the production of enough gametes [30]. 
It is therefore plausible that variation in reproductive success 
may be better characterized as a positive quadratic function of 
mating success. This is important because it suggests that sex 
differences may occur not only in the slope of linear functions 
of reproductive success (mating success), but also in the shape 
of such functions. The comparison between the maximum of 
the curvilinear function with the average mating success of a 
population enables one to test whether the reproductive success 
of members of one sex in a population is limited by the avail- 
ability of mating opportunities (maximum > average mating 
success) or by intrinsic reproductive costs (maximum < average 
mating success [30]). 

Third, polyandry makes the relationship between mating and 
reproductive successes more complex [9,31 -33]. With polyandry, 
the reproductive success of an individual male is the product of 
his mating success, the fecundity of his mates (mate fecundity) 
and the proportion of ova fertilized within the clutch of each 
female he mated with (paternity share). Because Bateman gradi- 
ents are measured as the slope of the univariate least-square 
regression of reproductive success over mating success, they 



capture both direct selection on mating success, but also indirect 
selection on other components of reproductive success if mating 
success covaries with mate fecundity and /or paternity share. 
Hence, the Bateman gradient behaves as a selection differential 
rather than act as a selection gradient as it encompasses direct 
and indirect selection on mating success [34]. For example, the 
male Bateman gradient can be inflated via paternity share if, for 
a given mate fecundity, males that mate with more females 
have on average a higher paternity share of each clutch [35] . Alter- 
natively, by mating with more females, males might suffer from 
reduced paternity share through trade-off mechanisms [36]. 
One way to measure the relationship between mating and repro- 
ductive successes controlling for these covariances is to use a 
traditional multivariate approach such as 

Reproductive success = /3X + 8, (1.1) 

where /3 is the vector of partial regression coefficients on mating 
success, mate fecundity and paternity share, X the matrix of 
phenotypic values for those traits and s an error term of 0 
mean [34]. 

Finally, it is becoming clear that females can also have 
positive Bateman gradients, even in species with typical 
'sex roles' [18,37-39], which suggests that sexual selection 
for mating success might also be significant in females 
[2,18,40-43] potentially owing to cumulative benefits associ- 
ated with mating [44,45]. Note that this is not necessarily 
inconsistent with the Darwin -Bateman paradigm provided 
that the Bateman gradient is steeper in males than in females 
in a given species. However, the significance of female Bate- 
man gradients and the relative strength of sexual selection on 
female mating success remain debated [22,45,46]. Alterna- 
tively, a positive female Bateman gradient may arise if 
more fecund females are exposed to higher mating rates 
solely because they are more attractive to males, thus rever- 
sing the causality of the relationship between female mating 
and reproductive successes [22,47,48]. Explicit tests of the 
significance of female Bateman gradients remain scarce. 

Here, we combine an experimental approach with multi- 
variate analyses to address these issues and resolve the 
significance of Bateman's principles in replicate groups of red 
junglefowl, Gallus gallus. First, we test Bateman's principles 
using the traditional approach of inferring mating success 
from genetic parentage of offspring (MSgen) and univariate 
regressions of reproductive success on mating success inferred 
from offspring parentage. Second, we measure mating success 
using fine-grained mating behaviour data to measure the bias 
introduced by the traditional use of mating success inferred 
from parentage. Third, we use a multivariate approach to quan- 
tify the relationship between male mating success and 
reproductive success, controlling for mate fecundity, paternity 
share and their covariances with mating success. Finally, we 
explore the causality of female Bateman gradients by testing 
whether variation in female mating success causes changes in 
female reproductive success. 

2. Material and methods 

(a) Observations and parentage assignment in 
semi-natural conditions 

We studied a population of red junglefowl at the field station of 
the University of Oxford between May -September 2007 and 



August -September 2008, as detailed in reference [35], Briefly, we 
observed 13 replicate groups of three adult males and four adult 
females in outdoor pens, typical group size and sex ratio in natu- 
ral groups [49]. In each replicate group, we monitored mating 
behaviour for 10 consecutive days. We collected all the eggs 
laid in each group from the second day of observation (i.e. the 
first day in which inseminations may have resulted in fertiliza- 
tion) and for the subsequent 10 days (i.e. day 2-11 inclusive, 
see the electronic supplementary material, table SI on variation 
in female reproductive success). One group departed from this 
pattern as egg collection occurred from day 7 to day 16. We 
tested the effect of removing this group from the analysis, but as 
no qualitative difference was observed, we present only the results 
of the fuU dataset. We incubated eggs artificially for 7 days when 
we collected the embryos of fertilized eggs. In order to measure 
whether sperm depletion was likely to create sexual selection in 
females, we recorded the percentage of eggs that showed no evi- 
dence of embryo development (indicating that the eggs were 
either not fertilized or suffered embryo death within the first 
hours of development) in five of the 13 replicate groups. In these 
five groups, 11% of the eggs showed no sign of embryo develop- 
ment. More than half of these eggs were laid in the first days of 
egg collection, suggesting that they were not fertilized owing to 
the expected delay between insemination and fertilization [50]. 
Excluding eggs laid in the first 2 days, reduced the proportion of 
eggs with no sign of embryo development to 5%. Collectively, 
these results indicate that if some fertilizations went undetected 
in the study, then these represented a small proportion of the 
eggs. All embryos (n = 254) were genotyped at seven variable 
microsatellite loci [50], their paternity and maternity were success- 
fully assigned in Cervus 3.0 ([35,51,52] and the electronic 
supplementary material, table S2). 

(b) General analytical approach 

Reproductive success in males and females was calculated as the 
number of embryos produced by each individual. Mating success 
inferred from genetic parentage of the offspring (MSgen) was 
calculated as the total number of mates with whom an individual 
produced offspring. Importantly, we also measured 'total' male 
and female mating success (MStot) by adding to MSgen any 
additional mate with whom a focal individual had been observed 
successfully copulating during the 10 day period of observation 
but to which no offspring were assigned. To control for the poten- 
tial effects of mating repeatedly with the same mate, we entered as 
a covariate the average number of successful matings that a focal 
individual was observed to have with each mate. As no difference 
was observed by adding this covariate, we did not report results of 
analyses that included it. In addition, for males, mate fecundity 
was calculated as the average number of eggs laid by all females 
with whom a male produced offspring, corrected by the number 
of eggs laid by all females with whom a male was observed to 
copulate successfully (see calculation of MStot above). Accordingly, 
the paternity share of a male was calculated as the proportion of 
embryos sired out of the total number of embryos produced by 
all the females with whom he successfully mated. We calculated 
the opportunity for total sexual selection I {Ij sensu [35]) as 
standardized variance in reproductive success [4,53,54]: 



7 — 2 

where cr represents the variance in reproductive success, and T 
is the square mean of the reproductive success of members of 
one sex in a group. Similarly, we calculated the opportunity 
for sexual selection on mating success. Is (Im sensu [35]) as the 
standardized variance in mating success 

h=^- (2.2) 
MS 



We also used a novel method described by Moorad & Wade [55] to 
determine the proportion of I explained by the different components 
of male sexual fitness, namely mating success, mate fecundity and 
paternity share. In brief, we partitioned I into additive components 
(either only Is in the univariate analysis, or with the other com- 
ponents of male sexual fitness in the multivariate analysis), and 
significance was tested by normal bootstrapping [56] to estimate con- 
fidence intervals. Finally, to obtain results comparable with other 
organisms with various population and clutch sizes, we produced 
a standardized measure of Bateman gradients [2]. We standardized 
Bateman gradients as usually done for selection gradients by divid- 
ing reproductive success by the average population score (relative 
fitness), whereas each individual trait (mating success, mate fecund- 
ity and paternity share) was subtracted by its population mean and 
then divided by its standard deviation in the population, to obtain a 
mean of zero and a standard deviation of unity. 

In the models comparing different ways of calculating the Bate- 
man gradients, estimates of all gradients were deduced from a 
simple linear model of male or female reproductive success and 
different predictors [5]. To obtain the significance of fixed effects 
and overall fit of models, we performed generalized-mixed 
models (kner, library lme4 in R v. 2.15.2, R Core Team [57]), with 
replicate group as random effect with 13 levels. Some birds were 
used in more than one replicate group (four males used in two 
groups, four in three groups, one in four groups, nine females in 
two groups and seven females in three groups). Repeatability in 
the components of reproductive success (total reproductive success, 
mating success, mate fecundity and paternity share) across groups 
was consistently low in males [35]. For males used in multiple 
groups within the same reproductive season, relative male's repro- 
ductive performances in one group were poor predictors of their 
relative performance in the next group (reproductive success: adj. 

= —0.07; mating success: adj. = —0.07; mate fecundity: 
adj. R^= -0.08; paternity share: adj. = -0.08). Similarly, the 
relative reproductive and mating success of a female in a group 
did not predict her reproductive success and mating success in 
the next replicate group within the same year (reproductive success: 
adj. R^ = —0.01; mating success: adj. R^ = —0.05), suggesting that 
the reproductive performance of these birds was largely contingent 
on the d3mamics of different replicate groups rather than consist- 
ently determined by inherent properties of the individual or by 
seasonal patterns. Nevertheless, bird identity nested within group 
was fitted as an observation-level random effect to control for 
pseudo-replication and overdispersion. In all mixed models, the 
response was male or female reproductive success which con- 
formed to a Poisson distribution; hence, we used a log link 
function to test the fixed effects. The significance of each fixed 
effect was tested by a log-Hkelihood ratio test comparing models 
fitted by the Laplace approximation including or excluding the 
tested fixed effect. However, because the log-likelihood ratio test 
is considered non-conservative for fixed effects [58], and because 
the log link function does not directly test the linear relationship 
that the Bateman gradient aims to describe, we calculated confi- 
dence intervals for each estimate of slope with bootstrapping 
methods. We first fitted a linear-mixed model (Ime, library nkne, 
in R V. 2.15.2) on reproductive success with the fixed effects. We 
then shuffled the residuals given by this model, added them to 
the fitted values and refitted the model to these newly created data- 
set [59]. To be conservative, we repeated this bootstrapping method 
2000 times on mixed models including either (i) the bird identity, 
(ii) the group of birds and (iii) the bird identity nested within the 
group of birds as random effects. The results presented include 
0.95 confidence intervals given by these simulations. 

(i) The traditional measure of Bateman's principles 

We calculated the total opportunity for selection (J), and the 
opportunity for sexual selection (Is) for males and females as 
the standardized variance in mating success inferred from 



genetic parentage (MSgen). We also calculated confidence inter- 
vals on I indices and the percentage of I explained by Is [55]. 
Male and female Bateman gradients were calculated indepen- 
dently by using male and female MSgen as fixed effects in the 
models above. To compare male and female Bateman gradients, 
we merged male and female datasets and computed mixed 
models fitted with a Poisson distribution that used, without 
intercept (i) sex and MSgen as fixed effects and (ii) sex, MSgen 
and the interaction between both as fixed effects. The comparison 
between (i) and (ii) by a log-likelihood ratio test quantified the 
statistical significance of the difference between male and 
female Bateman gradients. We also tested for quadratic effects 
by adding MSg^^^ as a fixed effect and calculated confidence 
intervals for MSgen and MSg^^^ as previously described. 



(ii) The bias introduced by inferring mating success 
from offspring parentage 

To investigate the effect of matings not resulting in fertilization 
on Bateman's principles, we recalculated male and female I, Is 
and Bateman gradients as above, using MStot instead of MSgen- 
We then compared the general fit and the Akaike information 
criterion (AIC) of models based on MSgen versus MStot- 



(iii) The multivariate measure of Bateman's principles 

To test whether the opportunity for sexual selection (Is) and the 
Bateman gradients were artificially inflated by the covariances of 
mating success with other components of male reproductive suc- 
cess [35], we added male mate fecundity and paternity share to 
the model. We first calculated opportunity of selection on mate 
fecundity and paternity share as described in reference [55]. 
Second, we tested the effect of adding mate fecundity and pater- 
nity share in the calculation of the Bateman gradient by 
comparing the goodness of fit and the AIC of the generalized- 
mixed models. We also tested the effect on male reproductive 
success of the covariances between all components of male repro- 
ductive success by adding second-degree interactions to the 
model, but as this was not significant we do not report it here. 

(iv) The significance of female Bateman gradients 

We explored the hypothesis that a positive female Bateman gradi- 
ent arises because mating with additional males causes an 
increment in female fecundity (adaptive hypothesis). The adaptive 
hypothesis predicts a specific temporal pattern, whereby female 
reproductive rate (i.e. the probability to lay a fertile egg on a 
given day) increases over successive days for polyandrous females, 
as they accumulate additional partners, but remains constantly 
low in monandrous females. There are two alternative (null) 
hypotheses for a positive female Bateman gradient. First, inher- 
ently more fecund females mate with more males. Second, even 
if all females mate with the same number of males, then inherently 
more fecund females will, in principle, by producing more eggs, 
have a higher probability of producing offspring sired by more 
males. However, measuring mating success based on behavioural 
data should make this latter explanation less likely to apply to our 
study. Neither of these alternative scenarios predicts temporal 
changes in the reproductive rate of polyandrous versus monan- 
drous females over the course of a trial. We therefore analysed 
variation in the probability that a female would lay a fertile egg 
on a given day over successive days of the trial in relation to her 
cumulative mating success (female MStot =1/ 2 or 3), with a 
mixed model fitting a binomial distribution and, including trial 
day, MStot, and day x MStot interaction as fixed effects. 

We then conducted an experiment to test the causal relationship 
between female mating and reproductive successes. Monoandry and 
polyandry treatments were sequentially conducted in random order 
on 13 females. Females received controlled matings [60] in the 



afternoon on 3 consecutive days, either three times with the same 
male (monoandry) or with three different males (polyandry). Egg 
laying rate was monitored for each of the females for the 7 days fol- 
lowing the treatment [50]. Females were housed in groups, and eggs 
were assigned to females using orally administered coloured lipid 
dyes [60]. We compared the number of eggs laid by a female follow- 
ing the polyandry treatment with that of the same female following 
the monoandry treatment, with a paired f-test. As our sample was 
rather limited (n = 13), we conducted a power analysis to estimate 
the likelihood of finding an effect size similar to the one observed 
in the semi-natural groups. We ran a model simulating a Poisson dis- 
tribution with the same sample size and average number of eggs laid 
as in the controlled mating experiment. The power analysis showed 
that we had 79.2% chance of detecting the difference in fecundity 
between females mated to one and females mated to two males 
observed in the semi-natural groups. Note that this power analysis 
is conservative because it does not consider the within-female 
paired design. 



3. Results 

(a) The traditional measure of Bateman's principles 

The opportunity for selection was almost twice as high in 
males as in females (table 1). Similarly, the opportunity for 
sexual selection on mating success was approximately four 
times higher in males than in females, and was significantly 
positive in both sexes (table 1). The variance in mating success 
explained 56.9% of the variance in reproductive success in 
males, but only 24.1% in females (table 1). Finally, the male 
Bateman gradient was significantly steeper than the female 
gradient, as indicated by the comparison between a model, 
including the fixed effects sex and MSgen, and a model, includ- 
ing sex, MSgen and the interaction between both (x^ = 9.4, 
p = 0.002; figure la). However, the Bateman gradient was 
significantly positive in both males and females (figure la 
and table 1). Standardized measures of male and female gra- 
dients produced a qualitatively similar pattern (table 1). We 
detected a positive quadratic relationship between repro- 
ductive and mating successes in females (estimate MSgen = 
5.65 [3.60; 7.60], estimate MSg,^ = -1.25 [- 1.89; -0.63], 
figure 2b) but not in males (estimate MSgen = 1-84 [-0.65; 
4.32], estimate MSg,^ = 0.61 [- 0.01; 1.23]). The female 
quadratic function was maximized when females had 2.26 
partners (female MSgen = 2.26 [CI 0.95; 6.04]), which was 
very close to the average female mating rate in the population 
(average MSgen = 2.14, figure lb). 

(b) The bias introduced by inferring mating success 
from offspring parentage 

Overall, 29.4% of the pairs that were observed mating did not 
sire offspring together and were therefore missed by estimates 
of mating success based on genetic parentage assignment 
(i.e. male and female MSgen). When taking into account mat- 
ings that did not result in fertilization by correcting data on 
mating success based on parentage with behavioural obser- 
vations (MStot)/ the opportunity for sexual selection on 
mating success dropped in males (male Is = 0.34 [CI: 0.16; 
0.53]) and became null in females (female Is = 0.03 [CI: 
— 0.06; 0.09]; table 1). Similarly, when correcting data on 
mating success based on parentage with behavioural obser- 
vations, the variance in mating success explained 42.5% of 
the variance in male reproductive success and only 5.4% 
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Figure 1. (a) Male (solid squares and line) and female (empty diamonds, 
dashed line) Bateman gradients based on MSgen- (b) The quadratic relationship 
between female MSgen and reproductive success over a MSgen range of 0-3, 
vertical bars represent CI. The function shows a maximum for MSgen = 2. 
The vertical dashed line shows (approximately) mean female MSgen across 
groups, the shaded section around the line represents SE around the mean 
(2.13 + 0.16). 



of the variance in female reproductive success (table 1). 
Correcting data on mating success based on parentage with 
behavioural observations also greatly changed the models 
estimating male and female Bateman gradients (table 1). The 
Bateman gradient remained significantly positive in both 
males and females (table 1) after correcting data on mating 
success based on parentage with behavioural observations. 
However, both male and female gradients were reduced, 
indicating that the traditional use of mating success based on 
offspring parentage can lead to considerable errors in the esti- 
mates of both unstandardized and standardized Bateman 
gradients for males and females (table 1). 
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Figure 2. Average number of eggs produced by females mated three times 
with a single male (monoandry) or once with each of three males (polyandry). 
(Online version in colour.) 

(c) The multivariate measure of Bateman's principles 

Taking into account the covariances between male components 
of reproductive success explained more than 96% of the oppor- 
tunity for selection in males (J), resulting in a significant 
improvement on the explanatory power of male reproductive 
success and in a lower AIC than that of the univariate model 
based on MStot (table 1). This multivariate model identified 
three independent sources of selection on male reproductive suc- 
cess: mating success (i.e. Bateman gradient), mate fecundity and 
paternity share (table 1), and revealed that, by failing to control 
for covariances between these components of male reproductive 
success, the traditional univariate approach overestimated the 
male Bateman gradient by more than 153% of the unbiased selec- 
tion gradient (table 1). Again, the use of standardized gradients 
qualitatively agreed with this result (table 1). 

(d) The significance of female Bateman gradients 

Contrary to the predictions of the adaptive hypothesis and 
consistent with those of the null hypotheses, we found that 
the probability of laying an egg did not change with time 
(i.e. throughout the duration of a trial), but was consistently 
higher for females with higher mating success (MStot- 
xi = 6.79, p = 0.01, day: xi = 0.02, p = 0.88, day x MS^t: 
xi = 0.080, p = 0.78). 

In the controlled mating experiment, females did not pro- 
duce more eggs when they were mated to different males 
(polyandry) than when they were repeatedly mated to the 
same male (monoandry, t = —0.41, d.f. = 12,p = 0.69; figure 2). 

4, Discussion 

Bateman's principles are a cornerstone of modern sexual selec- 
tion theory, yet intense recent debate has called into question 
their measure and relevance in studying sex differences 
[18,21,61-64]. Our study contributes to this ongoing debate 
by demonstrating that the traditional use of MSgen can generate 
severely misleading estimates of Bateman's principles. 
Because matings can fail to result in fertilization, particularly 
in polyandrous species, the number of females with whom a 



male mates successfully can only be as great as the number 
of females with whom he sires offspring and typically lower 
than this (i.e. MStot > MSgen). This inequality means that infer- 
ring mating success from offspring parentage is likely to 
introduce a systematic bias, by overestimating the opportunity 
of sexual selection on males. Importantly, for a constant share 
in paternity and mate fecundity, the total reproductive success 
of a male is directly proportional to his mating success inferred 
from offspring parentage but is independent from residual 
mating success that fails to result in paternity. Failing to con- 
sider this latter component of mating success is therefore 
likely to overestimate the steepness of the male Bateman gradi- 
ent [17], and, because the probability of fertilizing at least one 
egg increases with clutch size, this bias is likely greater in 
species with smaller clutches. 

Another large source of error arises from the covariance of 
male mating success with the other constituents of male 
reproductive success, namely mate fecundity and paternity 
share, which are not considered in the traditional univariate 
approach to Bateman's principles. In red junglefowl, male 
mating success and paternity share may covary for different 
reasons. In our study population, males with high mating 
success also tend to enjoy high paternity share. This appears 
to be mediated by social dominance which allows a male to 
access more females and repeated mating opportunity with 
each female, which conveys a sperm competition advantage 
(high paternity share), likely through a replenishment of the 
male's sperm in the female's sperm storage tubules [35]. An 
alternative scenario may arise when males invest in either 
mating opportunities or fertilization efficiency through 
alternative mating tactics [65]. Under such conditions, 
males with high mating success may be poor sperm compe- 
titors and suffer low paternity share, whereas males with 
high paternity share may suffer low mating success, resulting 
in negative (rather than positive) covariance between mating 
success and paternity share. Several studies have shown that 
in domestic fowl, socially dominant males tend to produce 
ejaculates of lower sperm swimming velocity than the ejacu- 
lates produced by their subordinates [66-68]. Controlling for 
mating frequency and number of sperm inseminated, higher 
sperm velocity would give subordinate males an advantage 
in sperm competition [69]. A negative covariance of mating 
success and paternity share may also be more likely to 
occur in species where repeated matings are less relevant, 
such as for example external fertilizers. Covariance between 
male mating success and mate fecundity may also introduce 
significant bias, if for example males with high mating 
success mated preferentially with more fecund females. 

Measuring Bateman gradients taking into account male 
paternity share and mate fecundity is non-trivial and to the 
best of our knowledge, few studies have attempted this in 
the past. Fritzsche & Arnqvist [39] have used the covariance 
between a trait and the residual reproductive success unex- 
plained by the Bateman gradient to measure the overall 
selection on mate fecundity and paternity share. This 
approach, however, is simplistic: the potential covariances 
between mating success and other components of male 
sexual fitness make reality more complex. We really need to 
integrate mating success, mate fecundity and paternity 
share in a multivariate framework that takes into account 
their covariances in order to accurately quantify sexual selec- 
tion on mating success (i.e. the Bateman gradient), mate 
fecundity and paternity share. Few studies have attempted 



to quantify such covariances. A recent study of the hermaph- 
roditic snail, Physa acuta, found significant covariances 
between different components of male reproductive success 
[28]. Similarly, studies of passerine birds forming social pair 
bonds have quantified the covariance between within- and 
extra-pair components of male reproductive success, demon- 
strating that such covariance can represent a substantial 
source of variation in male reproductive success in some 
populations [70,71], but less so in others [48,72]. In part, 
differences between these studies are likely to reflect a combi- 
nation of biological factors, such as population size, patterns 
of variation in polyandry, mate availability and clutch size, 
and future studies should seek to resolve how Bateman prin- 
ciples can be modulated by these factors. Some discrepancies 
may also reflect methodological differences. For example, 
Pelissie et ah [28] attributed variance in remating rates with 
the same female to sexual selection on male mating success 
[28], although in red junglefowl mating repeatedly with the 
same female has been shown to have a direct influence on 
paternity share [35]. Similarly, the distinction between 
within- and extra-pair reproductive success is not readily 
applicable to species that lack social pair bonds. 

Our results also confirm the concerns expressed by pre- 
vious studies over the interpretation of female Bateman 
gradients [47,48]. Positive female Bateman gradients, while 
significantly lower than corresponding male Bateman gradi- 
ents, have been shown in a number of species, for example, 
bank voles [15], wild turkey [16] and Drosophila [14,18]. How- 
ever, the causality of the relationship between female mating 
and reproductive successes is often unclear. Using a tra- 
ditional approach, one would be led to conclude that red 
junglefowl females are selected to mate with multiple 
males, and optimize mating success at around two males. 
Because this value does not differ significantly from the aver- 
age female mating success observed in these groups, one 
would also be led to conclude that females are in control of 
mating rates [30]. However, these conclusions are misleading. 
Such a positive female Bateman gradient is likely to arise as a 
spurious consequence of more fecund females mating with 
more males, rather than as an adaptive consequence of 
female polyandry. That the adaptive hypothesis is unlikely 
to explain our observation of a positive female Bateman gra- 
dient is indicated by (i) the lack of a time lag, which would be 
required for polyandry to cause an increment in female 
fecundity, and (ii) the lack of any evidence that mating 
with additional males causes an increase in fecundity in the 
pair-wise controlled experiment. The null hypothesis, that 
females that are inherently more fecund attract more mates, 
is entirely consistent with our results and also more plausible 
given the biology of the study species. This pattern may arise 
through both female- and male-driven mechanisms. Female 
fowl tend to display a higher propensity to mate in periods 
when they ovulate than in periods when they do not, and 
males have been shown to preferentially target ovulating 
females [73], indicating that females that lay more eggs are 
likely to be more sexually promiscuous. Kokko et al. [62] 
emphasize the importance of establishing causality to under- 
stand Bateman's principles. Our study indicates that this is 
particularly important when it comes to the interpretation 
of female Bateman's gradients. 

Clearly, the approach that we adopted to measure male 
reproductive success also suffers from its own limitations. 
First, we may have missed a proportion of matings. Although 



we monitored the populations throughout the daily peaks in 
mating activity (early morning and late afternoon [74]), it is cer- 
tainly possible that some mating activity went unnoticed. If 
these matings resulted in fertilization, then they would have 
been accounted for by our measure of mating success which 
complements behavioural data with cases in which paternity 
data identifies mating events undetected by behavioural obser- 
vations. In actual fact, such cases were rare in our study (7.3% of 
the total number of mates was missed by behavioural data and 
detected through offspring parentage assignment). This does 
not eliminate the possibility that we may have missed mating 
events that failed to result in fertilization. However, undetected 
matings that failed to result in fertilization would mean that the 
discrepancy between MSgen and MStot is even greater than 
reported by our study, and thus the bias introduced by inferring 
mating success exclusively from offspring parentage is conse- 
quently also even greater than our study estimates. Second, 
the use of behavioural data requires some careful consideration 
and may need to be tailored to the specifics of different study 
organisms. In our study, we considered only sexually mature 
males and females during their breeding season, and counted 
only mating events that comprised the entire succession of be- 
havioural steps concluding with a successful cloacal contact 
[35]. A more careful approach may be required by studies of 
natural populations with limited possibilities for experimental 
control, less complete information on the reproductive status 
of different individuals or where matings observed represent a 
non-random subset of the mating occurring in the population 
[75]. In addition, although we incubated eggs artificially and 
sampled embryos early in their incubation (on day 7 of 21 day 
incubation period), we cannot entirely rule out the risk that 
our measure of reproductive success may have been influenced 
by embryo mortality at very early stage of development (i.e. 
within the first approx. 24 h from fertilization), which may be 
difficult to detect without molecular assays. However, the pro- 
portion of eggs that was deemed infertile in our study 
represented approximately 5% of aU the eggs collected. There- 
fore, if it occurred, undetected early embryo mortality would 
appear to have had only a modest contribution. 

Finally, our study shows that while efforts to reduce 
error can generate estimates of Bateman's principles that are 
quantitatively drastically different from estimates derived 
from traditional approaches, the patterns of sex-specific differ- 
ences observed are qualitatively entirely consistent with the 
three Bateman's principles. Namely, we found higher I and 
Js in males than in females and stronger sexual selection on 
male rather than on female mating success. Therefore, while 
our study suggests that previous estimates of Bateman's prin- 
ciples are likely biased, it also provides a robust confirmation 
that these principles are real and not entirely artefacts of meth- 
odological limitations. In conclusion, our experiment confirms 
the validity of Bateman's principles in a semi-natural popu- 
lation, but simultaneously reveals the necessity to rethink 
traditional approaches to study the evolutionary ecology of 
sex roles. 
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